Study of CP{N — 1) 6'- Vacua by Cluster-Simulation of SU{N) Quantum Spin Ladders 
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D-theory provides an alternative lattice regularization of the (1 + l)-d CP{N — 1) quantum field 
theory. In this formulation the continuous classical CP{N — 1) fields emerge from the dimensional 
reduction of discrete SU{N) quantum spins. In analogy to Haldane's conjecture, ladders consisting 
of an even number of transversely coupled spin chains lead to a CP{N — 1) model with vacuum 
angle 6^ = 0, while an odd number of chains yields 9 = tt. In contrast to Wilson's formulation of 
lattice field theory, in D-theory no sign problem arises ai 8 — n, and an efficient cluster algorithm 
is used to investigate the S-vacuum effects. At 6 = tt there is a first order phase transition with 
spontaneous breaking of charge conjugation symmetry for CP{N — 1) models with N > 2. 



PACS numbers: 05.50.+q,11.10.-z,75.10.-b,75.10.Jm 

CP{N — 1) models are interesting (1 + l)-d quantum 
field theories [1] which share a number of important fea- 
tures with (3 -I- l)-d QCD. These include asymptotic 
freedom, the dynamical generation of a mass-gap, and 
an instanton topological charge leading to non-trivial 9- y 
vacuum effects. Despite the fact that Nature has chosen 
= 0, it is an interesting challenge to understand the 
structure of 0-vacua, which in QCD may, for example, 
lead to 't Hooft's oblique confinement phases [2]. It has 
been conjectured that CP{N — I) models have a phase 
transition at 9 = tt. In the CP{1) — 0(3) case this 
phase transition is known to be second order with a van- 
ishing mass-gap [3-5]. For iV > 2, on the other hand, 
the transition is conjectured to be first order [6], which 
is consistent with large N analytic results [7]. For finite 
A'^, however, the investigation of these nonperturbative 
problems is highly nontrivial. In particular, in contrast 
to 0{N) models and other (1 -I- l)-d quantum field the- 
ories, CP{N — 1) models with N > 2 cannot be solved 
with the Bethe ansatz. This is because in the CP{N — 1) 
case an infinite set of classical symmetries is anomalous 
and can thus not be used to solve the quantum theory 
analytically. 

The numerical investigation of lattice CP{N—1) mod- 
els is also far from being straightforward, even at 9 = 0. 
Again, in contrast to 0{N) models which can be studied 
with the efficient Wolff cluster algorithm [8] , no efficient 
cluster algorithm is available for CP{N — 1) models [9]. 
There is even a no-go theorem that forbids the construc- 
tion of an efficient Wolff-type embedding algorithm for 
these models [10]. Still, at = a rather efficient multi- 
grid algorithm was developed in [11]. However, at 9 = tt 
the situation is much worse due to a very severe sign 
problem: the contributions from odd topological charge 
sectors almost completely cancel those from even charge 
sectors. This makes it exponentially hard to access large 
lattices which is necessary for reaching reliable conclu- 
sions about the phase structure. For this reason, previ- 
ous numerical studies of 0-vacua were limited to mod- 
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FIG. 1. Spin ladder geometry: the open circles belong to 
sublattice A, while the filled circles form sublattice B. 



erate volumes [12,13] or rely on additional assumptions 
[14]. In the CP(1) = 0(3) case, a Wolff-type meron- 
cluster algorithm allows efficient simulations at 9 = tt [5]. 
Unfortunately, due to the no-go theorem this algorithm 
cannot be extended to higher CP{N — 1) models. 

In this paper, we present a method that allows us, for 
the first time, to perform accurate unbiased numerical 
simulations of any CP{N — 1) model at 9 = tt. Our 
method is based on the D-theory formulation of field the- 
ory [15] in which continuous classical fields emerge dy- 
namically from the dimensional leduction of discrete vari- 
ables. D-theory provides an alternative lattice regulariza- 
tion of field theory which also yields the same universal 
continuum theory but is otherwise very different from 
Wilson's lattice field theory. In the case of CP{N — 1) 
models, the discrete variables are generalized quantum 
spins T° = ^X% which generate an SU{N) symmetry 
[Tx^T^] = i^xyfabcT^- Here fabc are the structure con- 
stants of SU{N). The spins are located on the sites x of 
a square lattice with spacing a of size Lx L' , with L ^ L' 
and with periodic boundary conditions. Hence, as shown 
in figure 1, we have a quantum spin ladder consisting of 
n = L' /a transversely coupled spin chains of length L. 
The x-direction of size L corresponds to the spatial di- 
mension of the target CP{N — 1) model, while the extra 
j/-dimension of finite extent L' will ultimately disappear 
via dimensional reduction. We consider nearest-neighbor 
couplings which are antiferromagnetic along the chains 
and ferromagnetic between different chains. Hence, the 
lattice naturally decomposes into two sublattices A and 



B with even and odd sites along the x-direction, respec- 
tively. Note that neighboring sites along the transverse 
y-direction belong to the same sublattice. The spins T" 
on sublattice A transform in the fundamental represen- 
tation {N} of SU{N), while the ones on sublattice B 
arc in the antifundamental representation {N} and are 
thus described by the conjugate generators — T^*. The 
quantum spin ladder Hamiltonian is given by 
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where J > 0, and 1 and 2 arc unit-vectors in the spatial 
X- and y-directions, respectively. Note that, in the SU{2) 
case, the {2} and {2} representations are unitarily equiv- 
alent. In particular, by a unitary transformation one can 
turn -T^* into T^. This is not possible for TV > 2. By 
construction the system has a global SU{N) symmetry, 
i.e. [H, T"] ~ 0, with the total spin given by 
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which satisfies the SU{N) algebra [T",r''] = ifabcT''. 

A priori it is not obvious that the (2 + l)-d quantum 
spin ladder provides a viable rcgularization of the (1-1-1)- 
d continuum CP{N — 1) field theory. As a necessary 
prerequisite, the quantum spin ladder does indeed have 
the global SU{N) symmetry of CP(A^ — 1) models. Using 
the coherent state technique of [16] one finds that, at zero 
temperature, the infinite system (with both L, L' -^ cxd) 
undergoes spontaneous symmetry breaking from SU{N) 
to U{N — 1). It should be noted that the choice of the 
SU{N) spin representations (in this case {A^} and {A^}) 
has an influence on the breaking pattern. Due to sponta- 
neous symmetry breaking, there are massless Goldstone 
bosons — in this case spin waves — described by fields 
in the coset space SU{N)/UiN - 1) = CP{N - 1). The 
CP{N — 1) fields P{x, y, t) are Hermitean N x N projec- 
tor matrices, i.e. P^ ~ P, with TrP = 1. Using chiral 
perturbation theory, the lowest-order terms in the Eu- 
clidean effective action for the spin waves are given by 

S[P] ^ f dt f dx f dy TT{p',dyPdyP 
Jo Jo Jo 

+ Ps[dxPd^P + \dtPdtP] - -Pd^PdtP}. (3) 

Here /3 = 1/T is the inverse temperature, ps and p^ are 
spin stiffness parameters for the x- and t/-direction, re- 
spectively, and c is the spin wave velocity. The action is 
invariant under global transformations P' = flPil' with 
ri S SU{N). The last term in the integrand of eq.(3) is 
purely imaginary and is related to the topological charge 



Q[P] = 1 f dt f dx TiiPd^PdtP], (4) 
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which is an integer in the second homotopy group 
Ii2[CP{N — 1)] = Z, and thus y-independent. Hence, 
the y-integration in the last term of eq.(3) can be per- 
formed trivially. This yields i9Q[P] where the vacuum 
angle is given hy 9 — L'n/a — rnr. Here a is the lat- 
tice spacing of the quantum spin ladder and L' /a = n is 
the number of transversely coupled spin chains. Hence, 
for even n the vacuum angle is trivial, while for odd n 
it is equivalent to 6 ^ tt. The same topological term is 
also generated when one uses a single antiferromagnetic 
spin chain with larger representations of SU{N). When 
one chooses a completely symmetric representation with 
a Young tableau containing n boxes in a single row on 
sublattice A and its anti- representation on sublattice _B, 
the resulting vacuum angle is again given hy 9 = nir. 

While the infinite (2 -I- l)-d system has massless Gold- 
stone bosons, the Coleman-Hohenberg-Mermin- Wagner 
theorem forbids the existence of massless excitations once 
the y-direction is compactified to a finite extent L'. As a 
consequence, the Goldstone bosons then pick up a non- 
perturbatively generated mass-gap m ~ 1/S^ and thus 
have a finite correlation length ^. Interestingly, for suffi- 
ciently many transversely coupled chains, the correlation 
length becomes exponentially large 



S, oc cxp{4:nL'ps/cN) > L' , 



(5) 



and the system undergoes dimensional reduction to the 
(1 + l)-d CP{N - 1) field theory with the action 

s[p] ^ f dt f dxTi {^[d^pd^p + \dtPdtP] 

Jo Jo 9 c 

- nPdxPdtP}. (6) 

The coupling constant of the dimensionally reduced the- 
ory is given by 1/g^ = L' ps/c. The mechanism of dimen- 
sional reduction is well-known from quantum antiferro- 
magnets [17,18] and occurs in all D-thcory models. It 
has also been discussed for SU{2) quantum spin ladders 
[19]. The dimensional reduction of ladder systems has al- 
ready been used in the D-theory construction of the 2-d 
0(3) model at non-zero chemical potential. The corre- 
sponding sign problem has been solved with an efficient 
meron-cluster algorithm [20]. In this paper, we extend 
the D-theory construction to CP{N — 1) models, which 
allows us to simulate them reliably aX 9 ~ Q and tt. 

One advantage of D-theory is that it allows us to con- 
struct quantum field theories using simple discrete de- 
grees of freedom instead of the usual continuum clas- 
sical fields. In particular, the partition function Z = 
Trexp(— /3i/) of the SU{N) quantum spin ladder can 
be written as a path integral using a basis of discrete 
SU{N) spin states q G {u,d, s, ...} on sublattice A and 



q G {u,d,s,...} on sublattice B. For SU{2) this corre- 
sponds to the usual t and J, spins. These can be simu- 
lated with the very efficient loop-cluster algorithm [21,22] 
which can even operate directly in continuous Euclidean 
time [23]. As mentioned in [24], this algorithm extends 
to SU{N) in a straightforward manner. 

We have used this cluster algorithm to investigate if 
CP{N — 1) models with N > 2 have a first order phase 
transition at 6 = tt where the charge conjugation symme- 
try C is spontaneously broken. A natural quantity that 
suggests itself as an order parameter for such a phase 
transition would be the topological charge Q[P] which 
indeed is C-odd, i.e. QfP] = Q[P*] ^ -Q[P]. Note 
that C is explicitly broken for 7^ 0, tt. At 6* = tt it is 
not broken explicitly because then the Boltzmann weight 
exp{ieQ[P]) = (-l)"?!^! is C-invariant. The topological 
charge itself is defined only in the framework of the target 
continuum theory. In the discrete spin system another or- 
der parameter, which is also C-odd, is more easily acces- 
sible. In the basis of quantum spin states q, q we define an 
order parameter Q[q,q] by counting the number of spin 
flips in a configuration. Q[q,q\ receives a contribution 1 
if a pair of nearest neighbor spins along the x-direction, 
QxQx+i^ flips to another state, q'^q' ,p at some moment 
in time. A spin flip from q^q^^i to q'^q' ., on the other 
hand, contributes —1 to Q[q,q]. In the quantum spin 
ladder, charge conjugation corresponds to replacing each 
spin state q^ by q,^_^_i (which is the conjugate of q^_^i)- 
Indeed, Q[q,q] changes sign under this operation while 
the action remains invariant. 

We have used the cluster algorithm to simulate SU{N) 
quantum spin ladders for N = 3,4, and 5, with n = 
L'/a e {2, 3, ..., 7}. The spatial size L/a has been varied 
between 20 and 300, and the inverse temperature /3J was 
chosen between 15 and 200. Our simulations conflrm the 
existence of a first order phase transition with sponta- 
neous C-breaking at 6* = tt for all N > 2. As expected, 
there is no phase transition a.t 9 — 0. Figure 2 shows 
Monte Carlo time histories of Q[q,q] for S'[/(4) spin lad- 
ders with n — 3 and 4 which correspond to CP(3) models 
a.t 9 = TT and 0, respectively. For n = 3 one observes two 
coexisting phases with spontaneous C-breaking, while for 
n = 4 there is only one phase which is C-symmetric. For 
n = 5 and thus 9 = tt one again finds a first order phase 
transition as in the n = 3 case. Figure 3 shows the prob- 
ability distribution of Q[q,q] for an SU{5) spin ladder 
with n = L'/a = 7, which corresponds to a CP(4) model 
at = TT. The double peak structure again implies a first 
order phase transition. A first order phase transition is 
also observed in the CP(2) case. 

At this point, we have shown that D-thcory indeed ver- 
ifies the conjectured 0- vacuum structure of CP{N — 1) 
models. In order to demonstrate that we approach the 
continuum limit of an asymptotically free field theory 
we have determined the correlation length ^(n) (defined 




FIG. 2. Monte Carlo time histories of the order parameter 
Q[q,q] for the CP(3) model at 6 = (n = 4, L = 180a, 
I3J = 50; and e = TV (n = 3, L = 200a, (3 J = 80). 




FIG. 3. Probability distribution of the order parameter 
Q[q,q] fortheCP{4) model at 6 = n (n = 7 , L = ma, f3J^6). 
Due to limited statistics the two peaks are sampled unevenly. 

using the second moment method [10]) in the CP{2) 
case as a function oi n — L'/a which controls the cou- 
pling l/g2 = L'ps/c. One obtains ^(2) = 4.82(4)a, 
\{A) = 17.6(2)a, and ^(6) = 61(2)a, which indeed shows 
the exponential increase of the correlation length charac- 
teristic for an asymptotically free theory. 

In D-theory CP{N — 1) models at 6* = can also be 
obtained by dimensional reduction of an SU{N) quantum 
ferromagnet with the Hamiltonian 
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In this case, using a symmetric SU{N) representation 
with a Young tableau containing n boxes, the symmetry 
again breaks spontaneously to U{N — 1) and the low- 
energy effective action for the spin waves takes the form 

S[P] = I dt fd^x Tiipsd^Pd^P-^ fdT PdtPdrP] 
Jo J '^ Jo 

^l3psJd^xTT[d^Pd^Pl (8) 
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FIG. 4. Monte Carlo data for the universal finite-size scal- 
ing function F{z) of the 2-d CP{2) model at6 = 0. The black 
points represent D-theory data from an SU{3) quantum ferro- 
magnet at (3 J = 5.5 and 6, while the open squares correspond 
to data obtained with the standard Wilson lattice field theory 
at l/g^ = 2.25 and 2.5. 

where x now represents a 2-d space-time coordinate. The 
Wess-Zumino term with the quantized prefactor n in- 
volves an interpolated field P{x,t,T) (r G [0,1]) with 
P(a;,i, 1) = P(x,t) and P(a;,t,0) = diag(l,0, ..., 0). The 
Wess-Zumino term vanishes after dimensional reduction 
and one obtains the action of the 2-d CP{N — 1) model 
at 9 = with 1/g^ = f3ps- In order to verify explicitly 
that this defines the CP{N — 1) model in the continuum 
limit, we now compare physical results obtained with the 
SU{N) quantum ferromagnet and the Wilson formula- 
tion using the standard lattice action 



S[P] 



J2 Ti-[^xP. 



X+fl\ 
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The continuum limit is reached as g — > 0. A convenient 
physical quantity characteristic for a given model is the 
universal finite-size scaling function F{z) — ^(2L)/^(L). 
Here £,{L) is the correlation length in a finite system of 
size L (again obtained from the second moment method), 
and z = S,{L)/L is a finite size scaling variable that mea- 
sures the size of the system in physical units. Figure 
4 shows Monte Carlo data for F{z) obtained both from 
D-theory and from the standard Wilson approach. Up 
to small scaling violations, the agreement of the two 
data sets confirms that the SU{N) quantum ferromag- 
net indeed provides a valid lattice regularization of 2-d 
CP{N — 1) models. Thanks to the cluster algorithm, the 
D-theory framework allows calculations that are much 
more accurate than the ones using Wilson's approach. 

To summarize, D-theory provides us with a powerful 
algorithmic tool, an efficient cluster algorithm that even 
works at 9 = TT. Despite great efforts, this remains impos- 
sible within Wilson's framework. In fact, a no-go theorem 
forbids the construction of an efficient Wolff-type embed- 
ding algorithm. In the D-theory framework this theorem 



is evaded due to the use of discrete quantum variables 
instead of continuous classical fields. Also in Wilson's 
SU{N) lattice gauge theory all attempts to construct ef- 
ficient cluster algorithms have failed. The present study 
of (1 -I- l)-d CP{N — 1) models raises hopes that the D- 
theory formulation of (3 + l)-d QCD [15] could also be 
investigated with efficient cluster algorithms. In this way, 
even exotic phenomena like oblique confinement [2] may 
eventually become accessible to numerical investigation. 
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